Extended variational cluster approximation for correlated systems 
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The variational cluster approximation (VCA) proposed by M. Potthoff et al. [Phys. Rev. Lett. 
91, 206402 (2003)] is extended to electron or spin systems with nonlocal interactions. By introducing 
more than one source field in the action and employing the Legendre transformation, we derive a 
generalized self-energy functional with stationary properties. Applying this functional to a proper 
reference system, we construct the extended VCA (EVCA). In the limit of continuous degrees of 
freedom for the reference system, EVCA can recover the cluster extension of the extended dynamical 
mean-field theory (EDMFT). For a system with correlated hopping, the EVCA recovers the cluster 
extension of the dynamical mean-field theory for correlated hopping. Using a discrete reference 
system composed of decoupled three-site single impurities, we test the theory for the extended 
Hubbard model. Quantitatively good results as compared with EDMFT are obtained. We also 
propose VCA (EVCA) based on clusters with periodic boundary conditions. It has the (extended) 
dynamical cluster approximation as the continuous limit. A number of related issues are discussed. 

PACS numbers: 71.27. +a, 71.10.Fd, 71.15-m 



I. INTRODUCTION 



In the past decade, the dynamical mean-field the- 
ory (DMFT) was developed as a powerful tech- 
nique for the study of correlated electron systemsi^ 
Much knowledge have been obtained on the strong- 
correlation effect in three-dimensional systems, such as 



the Mott-Hubbard metal-insulator transition, 



itiner- 



ant ferromagnetism,5i& and the colossal magnetic resis- 
tance in manganitesj2, etc. Recently there are DMFT 
studies combined with ab initio energy band calculations 
for real materials&S One of the many ways to derive the 
dynamical mean-field theory is to approximate the exact 
Luttinger-Ward functional by a local one which depends 
only on the local Green's functions^ The latter is then 
obtained from an effective single-impurity model with a 
self-consistently ascribed bath. This approximation be- 
comes exact in the limit of infinite spatial dimensions. 1 
For finite-dimensional systems, this approximation ig- 
nores the spatial fluctuation effects from the outset, and 
the nonlocal interactions in the model are treated only 
on the Hartree level. 

In order to take into account the spatial fluctuations 
beyond DMFT, besides the systematic 1/D expansion 
approachjiSiii various cluster algorithms, e.g., the cellu- 
lar DMFT (CDMFT); 12 ! 13 : 14 ! 15 and the dynamical clus- 
ter approximation (DCA)^& have been proposed. The 
former uses clusters with open boundary conditions, and 
the latter has periodic boundary conditions. These meth- 
ods not only produce a k-dependent self-energy, but also 
include the short-range part of the nonlocal interaction 
within the cluster. On the other hand, DMFT has also 
been extended to treat the nonlocal interaction terms 
in the Hamiltonian more satisfactorily. In this respect, 
there is the extended DMFT (EDMFT) for including the 
nonlocal density-density or spin-spin interactions^iiSiiS 



and the DMFT for correlated hopping^! (DMFTCH) 
for treating the correlated hopping terms. An approach 
based on Active impurity models has been proposed to 
effectively circumvent the possible noncausality problem 
in previous cluster algorithms^ All these developments 
have led to various dynamical mean-field algorithms and 
their cluster extensions, 23 which constitute one of the 
most active fields of the strongly correlated electron the- 
ories. 

Recently, another cluster theory, the variational clus- 
ter approximation (VCA), was proposed by Potthoff et 
ui^L This theory provides a rather general framework 
to formulate cluster approximations for a lattice fermion 
systems. In the VCA, different approximations are real- 
ized by using different reference systems. Such reference 
systems can have different degree of complexity, ranging 
from a two-site system 2 ^ to a cluster embedded in several 
auxiliary baths. In the limit of continuous bath degrees 
of freedom, the VCA recovers the formula of CDMFT in 
the simple basis. The DMFT is obtained from the VCA 
by using a single-impurity reference system with contin- 
uous degrees of freedom. 

The VCA is constructed from a self-energy functional 
which is obtained by doing a Legendre transformation on 
the Luttinger-Ward functional $ [G] The Luttinger- 
Ward functional plays an important role in the theories 
of correlated fermion systems. It was first constructed in 
1960s in the course of expressing the thermodynamical 
grand potential through perturbative skeleton diagram 
expansions^ Later, it was shown that this functional is 
very useful in constructing conserving and consistent ap- 
proximations for correlated many-body systems The 
derivation of the VCA crucially depends on the following 
properties of the Luttinger-Ward functional^ (i) The 
self-energy as a functional of the full Green's function is 
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given by the functional derivation of $ [G] , 



(ii) When evaluated at the physical value of the Green's 
functions, $ [G] is related to the equilibrium grand po- 
tential fl through the relation^ 

f3n = $ [G] + Tr In [-G] - Tr (EG) . (2) 

(iii) $ [G] depends only on the interaction part of the 
Hamiltonian, and does not contain explicit dependence 
on the hopping term. In the derivation of the VCA, (i) 
and (ii) guarantee that the grand potential as a func- 
tional of the self-energy SI [E] can be constructed by Leg- 
endre transformation. This functional is stationary at the 
physical values of the self-energy E. (iii) guarantees that 
any reference system sharing the same interacting part of 
the original H will have the same form of the functional 
$ [G] , and hence the same form of the functional S7 [E] 
(although their actual values depend on the hopping part 
and may well be different). The VCA is constructed by 
identifying the stationary point of SI [E] in a self-energy 
function space limited by the reference systems as the ap- 
proximate solution of E. One of the merits of the VCA 
is that by using a reference system with small decou- 
pled clusters and limiting the number of variational pa- 
rameters, thermodynamically consistent and systematic 
physical results for a correlated electron system can be 
obtained with much less numerical effort 

One open question is, in the spirit of the VCA, how 
one can take into account the nonlocal interaction term 
beyond the Hartree approximation. In the VCA, it is 
possible to include longer- and longer-range nonlocal in- 
teractions by increasing the size of the cluster. How- 
ever, for interactions beyond the cluster size, employing 
the Hartree approximation amounts to defining a refer- 
ence system with a different nonlocal interaction from 
the original Hamiltonian, and hence with different func- 
tional SI [E] . Thus the condition for establishing the VCA 
is not fulfilled exactly and further approximation beyond 
the VCA is introduced. In this respect, the EDMFT 
was devised to incorporate the effect of the nonlocal in- 
teraction. The problems with EDMFT are that first, 
calculations are numerically very demanding because of 
the introduction of additional continuous bosonic bath(s) 
in EDMFT. Usually the quantum Monte Carlo (QMC) 
technique has to be used to solve the impurity problem, 
and hence the study is limited to finite temperatures^ 
Second, the EDMFT is basically a single-impurity ap- 
proach and can produce only local self-energies. To 
obtain the k-dependent self-energy, one needs a clus- 
ter extension of EDMFT such as the EDMFT+CDMFT 
scheme proposed by Sun et al.^ or the extended version 
of the DCA for including nonlocal interactions^ Both 
of them are very heavy numerical tasks. Therefore, it 
is desirable to extend the VCA to handle the systems 
with nonlocal interactions. In the continuous limit, the 



EDMFT and EDMFT+CDMFT should be recovered if 
one takes a single-impurity reference system and cluster 
reference system, respectively. Hopefully the flexibility 
and efficiency of the VCA will enhance the study of cor- 
related electrons with nonlocal interactions. Within the 
same framework, this kind of extended the VCA (EVCA) 
can be constructed for a system with correlated hop- 
pings. For this system, the EVCA theory in the con- 
tinuous limit should recover the DMFTCH if one uses 
a single-impurity reference system, and become a clus- 
ter extension of DMFTCH if one uses a cluster reference 
system. 

The purpose of this paper is to present such an ex- 
tension of the VCA. The original Luttinger-Ward func- 
tional is a functional of the full one-particle Green's 
function. To extend the VCA, we first express the 
grand potential functional SI [G, n] by the generalized 
Luttinger-Ward functional $ [G, n] i— This is accom- 
plished through a Legendre transformation. 36 Here, be- 
sides the one-particle Green's function G, a certain two 
particle-Green's function n should be chosen according 
to the form of the nonlocal interactions. For the construc- 
tion of $ [G, n] , we use a derivation that is formally dif- 
ferent from that of Chitra et al.^> but is more suitable for 
the present purpose. This is embodied, e.g., in that the 
universality of $ [G, n] and SI [G, n] is easily understood 
from the derivation process. Recently such a nonpertur- 
bative construction has appeared for the Luttinger-Ward 
functional <& [G] i2* In our work a similar construction for 
the case of multiple variables is carried out. 

Starting from Sl(G,n), the generalized self-energy 
functional^* Qevca [E,T] is obtained by doing another 
Legendre transformation with respect to G and n. The 
stationary value of this functional equals the exact grand 
potential SI of the original system^ and it is reached at 
the exact value of the self-energies. In contrast to SI [E] 
which is only independent of hopping term of the Hamil- 
tonian, the form of SIevca [E,T] is independent of not 
only the hopping term, but also the nonlocal interactions. 
It is hence possible to evaluate it in a clusterized refer- 
ence system sharing local interactions with the original 
Hamiltonian. Finally, the EVCA is obtained by identi- 
fying the stationary point in the generalized self-energy 
space of a reference system as the approximate solution 
of the original system. This part of our construction is 
in the same spirit as the VCAjSi although the derivation 
is formally different. 

In this paper, as specific examples, we describe the ex- 
tension of the VCA for two systems with nonlocal in- 
teractions. One is the Hubbard model with nonlocal 
density-density Coulomb repulsion, the other is the Hub- 
bard model with correlated hopping. In both cases, the 
Hamiltonians have a nonlocal interaction term that can 
be treated on the Hartree level within DMFT, and can be 
partly taken into account within the cluster in CDMFT 
or DCA methods. For the Hubbard model with density- 
density interaction, numerical results are given for the 
simplest reference system, in which each decoupled clus- 
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ter is a single impurity problem with three sites. We will 
call this scheme three-site EVCA. Its results are found to 
be quantitatively close to the EDMFT results. Finally, 
based on the real space formula of the DCA by Biroli et 
ai.|14 we propose a VCA and EVCA that utilize clusters 
with periodic boundary conditions. In the continuous 
limit, this version of the VCA will lead to the DCA in- 
stead of the CDMFT, and correspondingly, the EVCA 
formula will lead to the extended version of DCA. 

The plan of this paper is the following. In Sec. II, 
we present the construction of the EVCA for the Hub- 
bard model with nonlocal density-density interaction. In 
Sec. Ill, the EVCA is constructed for the Hubbard model 
with correlated hopping. Section IV is devoted to the 
EVCA formula with translation invariance. We discuss 
some issues related to the VCA and EVCA in Sec. V. A 
summary is given in Sec. VI. 



II. NONLOCAL DENSITY-DENSITY 
INTERACTION 



In this section, we will first obtain a grand poten- 
tial functional £1 [G, H] expressed via the generalized 
Luttinger-Ward functional $ [G, U] . The grand potential 
functional Qlw [G, IL] that bears stationary properties is 
obtained as a byproduct. The EVCA formula is estab- 
lished in general form. Numerical calculations are carried 
out for the three-site EVCA and results are tested. 

Our construction of $ [G, H] is nonperturbative and 
does not employ the diagrammatic expansion. Essen- 
tially, the way that we construct the multiple variable 
Luttinger-Ward functional is the same as that used by 
BaymSS for one variable case, but here a Legendre trans- 
formation is employed explicitly. Therefore, our deriva- 
tion relies neither on the applicability of the Wick's the- 
orem, nor on the convergence of the summation of the 
skeleton diagrams. It only relies on the assumption that 
the Legendre transformation, whenever required, can be 
carried out. This assumption has been argued to be cor- 
rect for the one- variable caseiS& 

Besides interacting one-particle Green's function, the 
generalized Luttinger-Ward functional is supposed to 
be also a functional of two-particle Green's functions. 
There are various two-particle Green's functions, such as 
the density-density Green's function, spin-spin Green's 
function, and the occupation-correlated electron Green's 
function, etc. Our way of constructing the generalized 
Luttinger-Ward functional is general and does not de- 
pend on the type of the Green's function. The type 
of the two-particle Green's function is determined by 
the specific form of the nonlocal interaction in the sys- 
tem. For the Hubbard model with a nonlocal density- 
density interaction, the generalized Luttinger-Ward func- 
tional <& is considered as a functional of the full Green's 
function G and the full density-density Green's function 



% (r 



(T T m (t) rij (t')). For other systems, it is 



particle Green's functions as variables of $. The Hub- 
bard model with correlated hopping will be considered 
in Sec.III. 



A. Generalized Luttinger-Ward functional 

We start from the Hubbard model with nonlocal 
Coulomb repulsion, 

i,j.cr i,j i 

+H loc . (3) 

Here, the notations arc standard. YV . indicates the sum 

— >j 

of sites i and j independently [i is the chemical potential, 
and Hi oc = !7 JV n^riii is the local Hubbard interaction. 
This model was recently studied using EDMFT: 34 The 
statistical action for Eq. (3) reads 

S[c*,c, G Q ,U ] 
= [ dr [ iT'^iWl-Go-Vfr-^^M 

J J ° ija 

+ f dr f dr' ]T m (r) [-I^ 1 ]^ (r - r') n, (t>) 

JO JO y 

+ / Hi oc {r)dT, (4) 
Jo 



where 



S(t-t>) (5) 



and 



n0±(T-T') = V ij 6(T-T') 



(6) 



The grand partition function 5 is expressed through 
the path integral as 



5 [Go, H ] = [ HVclVc^e- 3 ^^ ^, (7) 
J icr 

and the grand potential of the system is 



fi[G ,n ] = --ln5[G ,n ] 

p 



(8) 



straightforward to incorporate the specific type of two- 



Equation (8) in fact gives out the grand potential Vt as 
a functional of Go and Ho. In the following, we will do 
the Legendre transformation on the functional [Go, Hq] ■ 
However, since we have explicitly assigned values to Go 
and Ho in Eqs. (5) and (6), we introduce Go and Hq 
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as variables of the grand potential functional. For the 
clarity of the derivation, it is useful to distinguish a vari- 
ational variable and the physical value of it. In the fol- 
lowing, we use the symbol with a tilde on top to denote 
a quantity that can be varied, and the symbol without a 
tilde to denote the fixed physical value of that quantity. 
(l as a functional of Go and IIo is then defined as 



Go, IIo 



In 5 



Go, n 



(9) 



Similar to Eqs. (7) and (4), we define S and S as 

a [G ,n ] - J Y[VclV Cl „e-~ s l c *' c < do ^ (10) 



and 



S[c*,c, Go, IIo] 
r/3 rP 



= I dT I dT 'H c t M l~ G o %* (r r') c ja (/) 
Jo Jo lja 

+ / dr f dr'^'-ni (r) : [-^0% (r - r') : n 3 (r') : 
Jo Jo ij 

+ f H' loc {r)dr. (11) 
Jo 

In the above equation, we have singled out the Hartree 
contribution from the interaction and put it into the local 
part H' loc . We have : rij (r) := n, (r) — (rii) and 

H'loc (r) = H loc (r) + ]T p t [2n t (r) - (m)} , (12) 

i 

where pi — J2j Iq dr'[— ITq \j (r — r') (rij). Note that 
the [— ITq 1 ] in pi is used here as a fixed value Eq. (6). By 
definition, Q equals the physical value f2 when Go = Go 
and ITo = n . From Eq. (11), it is important to observe 
that the functional form of (l[Go, IIo] is determined by the 
form of H' loc as well as the specific way of introducing Go 
and ITo m Eq- (H)- It does not depend on the values of 
Go and IIo. 

We define the interacting one-particle Green's function 
G and the connected two-particle density-density Green's 
function II as 

Gija (T - T>) 

= -{T T c ia (T)c\(T'))s 



Vc* a Vc ia c ia (r) c* ja (r') 



xe 



-S c*,c,Go,n 



(13) 



% (r - r') 



(T T [n t (r) - (m)] [nj (r') - <«,•>]>§ 

, J-J-0 J irr 



G 



x [m (t) - (m)} [rij (t') - (rij)} e 



-S c*,c,G ,n 



.14) 



Note that G (r — r') is anti-periodic in (r — r') and 
II (r — t') is periodic, with the period /?. The Fourier 
transformation of them corresponds to functions of 
fermionic and bosonic Matsubara frequencies, respec- 
tively. With the above definitions, it is easily seen that 



G ija (t - r') = -j3- 



5n 


Go, ITo 




5 


Go\ 


(r' 


-r) 



(15) 



and 



-/?- 



% (r - r') 



In the following, we will abbreviate them as 



sn 


Go, n 




6 






-r) 



(16) 



and 



G 



n 



Go, n 



Go, IIo 



59. 



= -f3- 



Go, Ilf 



SGo 1 



-a- 



Go, IIq 



SUo 1 



(17) 



(18) 



As functionals of Go and IIo, G and II will get their 
respective physical values when Go = Go and IIo = IIo. 
Note that for nonzero Gg 1 and rig 1 , even in the case 
Hioc — 0, the calculation of G and II is a complicated 
many-particle problem and their values arc in general 
not equal to Go and IIo. 

To do the double Lcgcndrc transformation for 
f2[Go,IIo], we introduce an intermediate functional 

F[G,n], 



G,n 



= n 



Go, n 
sn 



-Tr 



- Tr 

Go, ITq 



= Q 



Go, IIq 



1. 



5Q Go, IIo 

~ i Gn 

SGn 1 



■IK 



+ -Tr 



(GGo 1 ) 



(19) 



5 



Here, the trace TrA — ^2 k n a A^a (iw n ) e lw " 0+ is carried 
out over space, time, and spin coordinates. The deriva- 
tive of F is obtained as 



SF G,n 



and 



SG 



SF G,n 



Gn 



m. 



G,n 



G,n 



(20) 



(21) 



On the right hand of Eqs. (20) and (21), we have ex- 
plicitly written Go and IIo as functionals of G and II 
which are regarded as independent variables. That this 
can be done is implied by our original assumption. From 
Eq. (19), the functional dependence of f2 on the one- and 
two-particle Green's functions is obtained as 



G,n 



G,n 



^Tr (GG 



i) 



The next step to construct a generalized Luttinger- 
Ward functional is to introduce the one-particle self- 
energy E and the quantity T. In the following, we will 
define them as 



There are several important properties of this func- 
tional, parallel to the properties of the original Luttinger- 
Ward functional $ [G] . (i) The functional derivative of 
$[G, II] with respect to G and to II gives E[G, n] and 
T[G, II], respectively, 



6$ G,n 



SG 



s$ G,n 



su 



= r 



(28) 



(29) 



This is easily obtained from Eqs. (25)— (27). E and T 
acquire the physical values E and T if they are evaluated 
at physical Green's functions {G,II}. Here, the physical 
meaning of T = ITq 1 + oTI^ 1 is somewhat obscure. In 
analogy to the one-particle self-energy, it can be regarded 
as an effective density cumulant that plays the role of the 
self-energy for two-particle Green's functions^ (ii) The 
grand potential functional can be expressed by $ as 



G,n 



= <p 



G,n 



TrlnG-Tr ( GG^ 1 



aTrlnn-Tr (nn^ 1 



(30) 



r 



Go 1 



no 1 



G" 



oTT 



(23) 



(24) 



Here, the inverse should be considered as the matrix in- 
verse in the space, time, and spin coordinates. In the 
definition of T, we introduce a real parameter a. Phys- 
ically, the bosonlike nature of the operator may ad- 
vocate the value a = 1/2, and it was actually taken by 
many authors. 33 However, since rii (r) is not a strict real 
boson field, the value of a cannot be determined to be 
1/2 from first principle. We will discuss this issue fur- 
ther in Sec. V B. Combined with Eqs. (20) and (21), the 
above definitions lead to the following equations: 



S[F [G,n]/3-TrlnG] 
SG 



S[F[G,IL]/3 + aTr Inn] 
SU 



= r 



(25) 



(26) 



With Eqs. (23)— (26), we are now ready to intro- 
duce $[G,n], which is a generalization of the original 
Luttinger-Ward functional $[G]. It is defined as 



4> 



G,n 



(3F 



G,n 



Tr In G + aTr In n. 



(27) 



Here, G and n o are both considered as functionals 
of G and H Eq. (30) is obtained by combining Eqs. 
(22) and (27). Unlike the original Luttinger-Ward func- 
tional <& [G], <&[G, n] is nonzero in the case of Hi oc = 0. 
Therefore the explicit expression of f2[G, n] in this limit 
is still a nontrivial many-body problem, (iii) The form of 
the functional <&[G,n] is determined only by the form of 
Hi oc once the type of the two-particle Green's function n 
is chosen. This is easily understood as a consequence of 
our derivation process. As stated above, the functional 
dependence of f2[G ,n ] on {Gg 1 ,^ 1 } [Eqs. (9) — (11)] 
is determined only by the form of Hi oc as well as by the 
specific form of Eq. (11) in which Gq 1 and IIq 1 are intro- 
duced .^This^is alsq^true for the functional derivatives of 
fi[Go,no], G and n, and their inverse functionals. Fol- 
lowing Eq. (19), i^[G,n], the Legendre transformation 
of ^[Gq^IIq 1 ] also bears the same property. And our 
conclusion follows Eq. (27) immediately. 

In the original paper of Luttinger and Ward, a func- 
tional expression of fl that is stationary at the physical 
value of G is obtained. It is very important for devising 
various kinds of thermodynamically consistent approxi- 
mations. Here, the functional form of the grand poten- 
tial functional that is stationary at the physical value of 
{G, n} is obtained from a modification on f2[G, n] of Eq. 
(30); namely, in Eq. (30), we replace the Go and Ho with 
their respective physical values Go and Ho defined in Eqs. 
(5) and (6). We denote the resulting functional as Qlw- 
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It reads 



00, 



LW 



G,n 



= $ 



G,n 



TrlnG-Tr (GGn 1 



- aTrlnfi- Tr UUn 1 



(31) 



The functional Q,lw[G, II] satisfies the stationary prop- 
erties 



50,lw g, n 



SG 



iG=G,n=n 



5Q 



LW 



G,n 



STl 



iG=G.n=n 



= 0, 



g = g, n = n =n 



(32) 



The above derivation, when reduced to the case of one 
variable G, will lead to the functional Olvk[G]. It is 
actually what Luttinger and Ward obtained by skeleton 
diagram expansion, and can be taken as a starting point 
of constructing consistent approximations. However, for 
the construction of the EVCA in the following, we will 
not make use of the form of f2 lw [G, II] . Instead, we take 
the functional form fi[G, II] as a further starting point 
toward EVCA. 



B. EVCA for density-density interaction 

In this section, we continue to derive the EVCA for 
the Hubbard model with nonlocal density-density inter- 
action (3). The derivation is based on the generalized 
Luttinger- Ward functional <&[G, n] obtained in Sec. II 
A. The grand potention has been expressed by <&. To 
go further, we follow the original idea of Potthoff^i and 
apply a Legendre transformation to $[G, II], 



.4 



E,r 



G,n 

5$ 



-Tr 



Tr 



G,n 



G,n 



SG 



G 



G,n 



STl 



Tr 



n 



(£3) 



Tr (rn) . (33) 



For the intermediate functional Af£,rl, we have 



5A 



E,r 



ST, 



5A 



5T 



= -G 



n 



E,r 



s,r 



(34) 



With A[E,r], the Q as a functional of E and T is 
obtained £1 



E,r 



.4 



E,r 



Tr In G, 



oTr In 



(no 1 



(35) 



Here, {G 1 ,I1 1 } is considered as the functional of 
{E,f}. 

To construct the EVCA, however, we need a grand 
potential functional of E and T which reaches its physical 
value at a stationary point. The above functional J7[E, T] 
doesnot satisfy this requirement. Therefore, similar to 
the construction of f2z,iy[G,n], W e suggest to use the 
following functional to establish the EVCA, 



I3&EVCA 



A 



E,T 



E,T 



Tr In G, 



oTr In 



(no- 1 



(36) 



In Eq. (36), instead of using the variable {Go, Ho} as in 
Eq. (35), their physical values Eqs. (5) and (6) deter- 
mined by the original Hamiltonian are used. It is easy 
to verify that this functional has the desirable stationary 
properties, i.e., 



5Q E vca 



e,t 



ST, 



sn 



EVCA 



E,r 



ST 



is=s,r=r 



i£=s,r=r 



0, 



and2£ 



^evca 



E = E, r = r 



= n. 



(37) 



(38) 



Another important merit of 51 evca is that its functional 
dependence on {E, T} is only determined by the form of 
Hioc, and it has nothing to do with the concrete values 
of E and T, or equivalently of Go and Ho. This is the 
consequence of the Legendre transformation applied to 
$[G, n] which also has such a property. Combining Eqs. 
(35) and (36) to eliminate A, we can express Qevca by 



EVCA 



n 



E, T 
1 



Trln^-E) -a-Trln^ 1 



1, 

0' 



-Tr In 



(Go" 1 



a^Trln (h" 1 



- T 
(39) 



This equation plays a fundamental role in the construc- 
tion of EVCA. Based on the same functional for the ref- 
erence system, we can introduce the EVCA. 
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To conveniently express the EVCA equations, we first 
define the reference system, apply Eq. (39) to the ref- 
erence system, and finally introduce the approximation 
that leads to the EVCA equations. The reference sys- 
tem is composed of identical clusters of original lattice 
sites, and each site in the cluster may be coupled to a 
number of noninteracting bath sites. Hoppings and non- 
local interactions are allowed only inside each cluster, 
while the local interaction Hi oc on each site is the same 
as in the original Hamiltonian. For clarity, we introduce 
a subscript c to denote a matrix with zero matrix ele- 
ment^between sites on different clusters. For an exam- 
ple, G C ij = if sites i and j belong to different clusters. 
Therefore, we can always write matrix G c in block diago- 
nal form by numbering the lattice sites in an appropriate 
order. denotes the submatrix on the Ith cluster. To 
project out the block diagonal part of a general matrix 
A, we use the symbol \A\. That is, \A\ij = Aij if i and 
j belong to same cluster, and \A\ij = otherwise. The 
submatrix of A on the Jth cluster is denoted by \A\ ! . 

The action of a general reference system is written as 

S ref [c\c,G 0c ,n 0c } = J2Sc[c^c,Gl c ,n I 0c } (40) 



The Green's functions are obtained as 



and 



s c [c*, c ,(4,i4] 

rl3 r & 



j drj dT'J2^A^[-Goc%Ar-r')c j(T (r') 



a r 

dr / 
I) Jo 



-'E 

ijei 



rfr'>-:n 4 :(r) [-I^ c % (r - r') : n 3 : (r') 



+ 



/ H'llc ( T ) dr. 
Jo 



(41) 



Here, H' t * c = UJ2 ie i n n n il + i 2n i 0) - ("»>]. : 

m := rii-irii), andp 4 = J2 3 J ^'[-H-t^h ( T ( n j)- 
Here and in the following, we use S c , H c , and tl c to de- 
note the corresponding functionals on each cluster. Since 
the clusters are geometrically identical, the functional is 
of the same form on each cluster, although the values of 
their variables could be / dependent. The partition func- 
tion and the grand potential for the reference system arc 
then written as 



Goc, n. 



0c 



i 



I,iEl,a 



(42) 



and 



Gqc, IIqc 



(43) 



G c 



Gqc, IIqc 



= -/3- 



Sil GqcHqc 



and 



SG 0c 



Sfl GqcIIoc 



(44) 



(45) 



The previous derivations are done for arbitrary func- 
tions Gg" 1 and Ilg . They also hold for the reference 
system where Gq 1 = G^ 1 and IIq 1 = H^ 1 . In particu- 
lar, the forms of the functionals are not dependent on the 
variables. Therefore, one can use the same set of sym- 
bols of the functionals and simply replace the variables 
G, S, etc., with G c , S c , etc., in each previous expression. 
When expressed in the domain of the self-energies of the 
reference system, the functional Qevca[E,F in Eq. (39) 
becomes 



&EVCA 

= n 



-Trln (G^-E, 



Trln (G 1 - T c ) - ^Trln (Eg" 1 - £ c ) 



+ ^Tr In 



(V 



(46) 



Since the degrees of freedom in the reference system are 
decoupled between different clusters, the functionals can 
be expressed as a summation over clusters. For example, 
we have 



n 



G c ,n c ] =]To c [clni 



(47) 



and 



5, 



G ~ , IT 



TrlnG^-Tr (g^Y) 



-aTrln^-Tr (n^Y) . 



(48) 



Here the trace is carried out in the cluster subspace. To 
be complete, we write down the intermediate expressions 
for the reference system. 



6<S> C 



SGI 

g£, n£ 



SUl 



f^ = [n -Y + a[n-Y. (49) 
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Similarly, the generalized self-energy functional 
f2[E c ,r c ] reads, 

pn [e c ,F c 
i 

= 2[s c ,f c ] -^Trln([G - c Y-^) 
i 

+a^Trln([n- c Y-r^). (50) 



The functional A[E C , r c ] is also a summation, 



A 



and A r satisfies 



s C) r c ] = J> C [eJ,i^' 



(51) 



yl r l 



5A C 



sn 



-Gi 



-Hi 



yl V l 
^C! 1 C 



y, 1 r J 



(52) 



The EVCA requires that the stationary properties 
of fWcAp,?], Eq. (37) and (38), be satisfied in 
the subspace of E c and II C , in which the functional 
QevcaI^c, T c ] can be exactly obtained by solving the 
reference system. At this point the strategy of the EVCA 
is the same as the VCA. The EVCA equation is then ob- 
tained by identifying the approximate self-energies of the 
original system as the ones at the stationary point of the 
functional Qevca[E c , r c ]. It means that we require 



SQ 



EVCA 



SE 

SQevca 



e c , r c 



and 37 



Qapp = QEVCA 



5j C — Hqdh ,r c — r a 



^app-) Tc - r ^ app 



(53) 



(54) 



where E app , T app , and fl app are the EVCA results for 
each quantity. For a reference system with continuous 
degrees of freedom, by introducing Eq. (46) and making 
use of Eqs. (47)— (52), the above equations reduce to 

Gi = 



G 1 - E, 



n 



n " 1 - r« 



r-li/ 



(55) 



Equation (55) holds for each cluster index I, and the 
Gi and IL[ are determined by [G^Y and [n^Y of the 
reference system via Eqs. (40) — (45). Since the matrices 
on both sides of Eq. (55) are block diagonal, the cluster 
index I can be dropped. 

This set of equations apply in the case where the refer- 
ence system has continuous degrees of freedom, such as 
a cluster where each site is coupled to continuous baths. 
It can be regarded as an extensions of the EDMFT to- 
ward the cluster algorithm,— Note that in this formula, 
the clusters of the reference system have open boundary 
conditions like CDMFT. We therefore denote this exten- 
sion as EDMFT+CDMFT. It will generally lead to self- 
energies that are not spatially translation invariant. As in 
CDMFT, the self-energies for the lattice model should be 
reevaluated from the cluster self-energies obtained from 
Eq. (55). Different schemes of estimators have been dis- 
cussed in Ref. 13. When the solutions are confined to be 
spatial translation invariant, each cluster becomes iden- 
tical. In particular, for the reference system composed of 
one impurity embedded in two baths described by local 
Weiss field Gq 1 and Ilg 1 , the original EDMFT equations 
will be recovered. In this case, the right hand of the first 
two lines of Eq. (55) can be evaluated in momentum 
space, and the local Green's functions are written as in- 
tegrals with free density of states. The EVCA equations 
reduce to 



G (iw n ) 



n (iu n ) 



Po(e) 



iu> n + \i — e — E (iu n ) 
Vo(e) 



-de, 



-de, 



-e - T (iv n ) 
£ (iu> n ) = G^ 1 (iw n ) - G^ 1 (iuj n ) , 
r c (iv n ) = IT 1 (iv n ) + all -1 (iu n ) . (56) 

Quantities here are the local ones, and po(e) — 
l/NJ2 k 6(e-e k ), V Q (e) = l/NJ2 k S(e-V k ) are the den- 
sity distributions of the eigenvalues of hopping matrix 
and density-density interaction matrix. u> n and v n are 
the fermionic and bosonic Matsubara frequencies, respec- 
tively. £fe and Vk are obtained from Fourier transforma- 
tion of the hopping and Coulomb interaction parameters 
in the original Hamiltonian, 



= E-** 
= E^< 



-ik(i*£— rj) 



(57) 



l ~3 



If the reference system has only finite degrees of free- 
dom, the variation of {E c ,r c } cannot be arbitrary. We 
can use a finite number of variational parameters t and 
V to specify the Green's functions Gq c (t) and TLq c (V), 
respectively. The stationary point of Qevca[E Ci T c ] may 
be unreachable in the space of {t, V}, which corresponds 
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to a certain subspace of {E c ,r c }. In such cases, we 
identify the approximate self-energies of the original sys- 
tem as the ones at the stationary point of the functional 

h EVC A[z c (t,v),f c {lv)}, 



sn 



EVCA 



(t,v) ,f c (t, v) 




t — i-app i — ^app 



t=t av „,v=v c 



app ) v — » app 



(58) 



It further simplifies into 



TV 



G c 



n c + a 



(Go 1 



l SE t 



(IV) 



St 

sf, 



(IV) 



St 



(59) 



TV 



G c 



G n 



<5S, 



(IV) 



+Tr 



Il r + a 



n - 1 - r t 



sv 

-\ sf c (t, v) 



5V 



= 0. 



(60) 



Equation (58) is only one of the many possible ways to 
construct the discrete version of the EVCA equations. 
In Eqs. (59)— (60), the fermion and boson contribu- 
tions are combined through the cross derivation terms 

ST C (t, V^j /St and SE C (t, V^j /SV. We can also conceive 

another form of discrete EVCA which has the same con- 
tinuous limit as Eq. (55), but without the cross derivative 
terms. For example, we propose the following discrete 
version of EVCA which separates the contributions from 
fermions and bosons. Define 



= fi 



-Trln (Go. 1 - E c ) - ±Trln (g^ 1 - E c ) ; 



0' 



s c ,r c 



E C) Ic 



jTrln (n^ 1 - T c ) + ^Trm (n^ 1 - r c ) , (61) 



(3 V uc V p- 
we get a new version of the EVCA equations 



SO. 



(t, v) ,r c (t,V 



St 



sn h 



e c (t,v) ,r c (t,v) 



sv 



\v=v app = 0. (62) 



This approximation scheme leads to the equations 



Tr 



TV 



G c 



n c + a 



G 1 - E« 



(no- 1 



-r, 



sz c [t, v 
si 

ST C (t, V 



SV 



= 0. 



(63) 



Both Eqs. (59) and (60) and Eq. (63) are extensions of 
the VCA24 to include the nonlocal density-density inter- 
actions beyond Hartree approximation. They recover the 
continuous version Eq. (55) in the continuous bath limit. 
It is noted that Eqs. (59) and (60) are formally exact, 
while Eq. (63) is not. Equation (63) can be regarded as 
a further approximation based on Eqs. (59) and (60) 
to neglect the cross-derivative terms. In our numeri- 
cal study of the three-site EVCA, however, it is found 
that due to the very large energy scale difference between 
fermion and boson contributions to the grand potential, 
in the first version of EVCA, the cross-derivative terms 
will spoil the self-consistent equations and make it dif- 
ficult to find the stationary point. The second version, 
Eqs. (61) — (63), works much better. As shown in the 
next section, it gives out rather satisfactory results for 
the extended Hubbard model. 

In the above theories, it is rather flexible to select the 
reference system according to the nature of the problem 
at hand. One can select the size and the shape of the 
cluster ) 31 ! 32 the number of discrete bath sites, as well 
as the definition of the variational parameters used in 
the reference system. For comparable complexity of the 
reference system, it is expected that the lower spatial di- 
mensions the system has, the more efficient it is to use 
larger size of the cluster instead of to use larger number 
of bath sites. In the VCA study of the Mott-Hubbard 
transition, 24,39 qualitatively correct results have already 
been obtained by using a reference system in which each 
decoupled cluster contains only two sites, one impurity 
site and one fermion bath site. This shows the amaz- 
ing efficiency of the VCA. It is expected that EVCA 
may share these advantages of the VCA, and becomes a 
powerful technique for the study of strongly correlated 
electron systems. Interesting physical problems, such 
as different phases in the extended Hubbard model, the 
competition between the Kondo and Rudcrman-Kittel- 
Kasuya-Yosida (RKKY) interaction, etc., may be studied 
by EVCA using small reference systems. 



C. Three-site EVCA: An example 

In this subsection, as an example, we present some nu- 
merical results for the simplest realization of the EVCA 
for the extended Hubbard model. Chitra et al. stud- 
ied the influence of the nonlocal Coulomb repulsion on 
the Mott metal-insulator transition using EDMFT and 
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the projection technique^ They found that the long- 
range Coulomb repulsion will change the second-order 
Mott transition at zero temperature into a first-order 
transition. Here we use the second version of the dis- 
crete EVCA, Eqs. (62) and (63), to illustrate the im- 
plementation of EVCA, and focus only on the properties 
of the metal and the insulator. For simplicity, we only 
consider the half-filling and paramagnetic phase. More 
detailed studies on the extended Hubbard model will be 
published elsewhere. The Hamiltonian that we will study 
is Eq. (3), 



(64) 



Gq 1 and 1 are therefore the same as in Eqs. (5) and 
(6) . We consider a reference system composed of discon- 
nected lattice sites, each of which is coupled with two 
bath sites: one fermionic site and one bosonic site. The 
unit cell of the reference system is schematically shown 
in Fig. 1. This is an extension of the two-site DMFT 
that works well on a qualitative level in the study of the 
Mott-Hubbard transition^ 




u, n 



FIG. 1: Structure of a three-site unit cell of the reference sys- 
tem Hamiltonian (66) . Filled dot denotes impurity site, open 
dot and square denote fermion and boson bath site, respec- 
tively. 



In the following part of this section, we will drop the 
tilde on the top of the variables G, n, etc., since this 
will not introduce confusion. For the reference system 
described above, the action Eq. (41) on a single site is 
written as 

a Jo Jo 

r/3 r/3 

+ dr dr'nir) [-Tig 1 (r - t')\ n(r') 
Jo Jo 

r/3 rP 
+ U drn^ (r)ni (r) + 2p / drnir) — /3p(n) . 
Jo Jo 

(65) 



Here p — [1/no (ifo) + Ylj Voj]{ n )- The Hamiltonian for 
the three-site unit cell of the reference system (Fig. 1) 



can be written as 

H imp = e V] dicier + V V] (ajcg + 4 a o-) + U n 1 n l 



(H - 2p) E 4c + &b + 7 (fit + b) 



X c - C - 



(66) 



Here n a = 4 c CT . c CT and a a are the annihilation operators 
on impurity and fermion bath sites, respectively, b is 
the boson annihilation operator on the boson bath site. 
Integrating out the bath degrees of freedom in Hi mp , and 
comparing the resulting effective impurity action with 
Eq. (65), we obtain the following relations: 



G 0o f (iu> n ) =iu> n +fi- 

n o («*) = 



v 



(67) 



CI 



ref 



~ In + | In (l + e"*) 

p p 



ln(l - e~^) -p{n). 



(68) 



The expression for the functional Clf and fib in Eq. 
(61) now reduces to^Si 



1 

N~ 



n 



f 



4 +°° r 

= ^W-^X J <^ e Po(e)ln|l + G(iw n ) [A(iu>„) - e]| , 



ref 



2a 

T 



+oo „ 

E J deV (e)ln 



1 Il(iv n ) [$(iv n ) - e] 



2 | de^ln 



l--n(twb)[$(*wb)-e] 



(69) 



where 



A(iu n ) 
$(ii/ n ) : 



sw n - e 



(70) 



The Green's functions G and n can be obtained 
through the Lehmann expression after the three-sites 
Hamiltonian Eq. (66) is diagonalized numerically. They 
are 



n /, \ 1 \ - I H4 
G CT (twj) = > — 



|(n|4|m)| 2 



( e -/J^m + e -^») 



(71) 



11 



and 



n (in) (i?o) 



|(»M TO )I 2 ( C -PE* 



- E n — E„ 



-/3E n 



(72) 



n (ivo) 



E 



(n\n\m)\ 2 e-^ 



mn,(E m =E rl ) 



1 \- \(n\n\m)\ 2 I -0E„ 

E -E ^ 



(73) 



Here ^^mp — Trc ^ 1 



The densities of states of the matrices hopping t%j and 
interaction Vy are selected to be of semicircular form, 



Po(e) = 
Vb(e) = 



7rJT 2 
2 



Vw 72 



i^ 2 



(74) 



stationary with respect to the fermion parameters e, V, 
and 1/NQb with respect to the boson parameters £,7. 
The Green's functions at such stationary points are our 
solutions for the three-site EVCA. For the boson opera- 
tor b, we take Nb < 00 optimized orthogonal states as the 
basest The final results should converge with respect to 
Nb. For the half-filled system that we consider here, the 
particle-hole symmetry condition of the reference Hamil- 
tonian Eq. (66) is 



U 2 7 2 n 
N b = 00. 



(75) 



It is found that with Nb = 8, the particle- hole symme- 
try is already satisfied very well, and the results do not 
change on further increasing Nb- With this Nb and con- 
sidering the conservative quantity n-f and nj_, the size of 
the largest matrix to be diagonalized is 48. This enable 
us to scan the parameter space quite quickly. For the 
parameter a that is unfixed within our theory, we take 
the value a = 1/2, same as in other EDMFT studies. 33 
The influence of different selection of a will be discussed 
elsewhere. 



These densities of states can be realized by nearest- 
neighbor couplings tij and on the Bethe lattice in the 
limit of the coordination number z — > 00, with scaling 
t — > W/2,/z and V — > Vq/2s/z, respectively. It is noted 
that within the present formula, the description of the 
nonlocal Coulomb repulsion Vij is only through two quan- 
tities V Q (e) = l/NY< k $ (Vfc - e) and £\ V oj = V k =o- It 
is apparent that the description is far from complete. A 
complete description should take into account different 
contributions from each k component of Vij, which can 
be realized partly by using larger a cluster in the refer- 
ence system. 

Another point is that using the above scaling and 
in the limit z — > 00, we get a diverging Vk=o- This 
reflects the fact that our present EVCA formula for 
density-density interaction is not a well-defined theory 
for infinite-dimensional models. The same problem exists 
for the symmetry-broken EDMFT. However, this doesnot 
prevent the theory from being a good approximation for 
finite-dimensional systems. In our calculation, we simply 
take Vfc = o as an independent parameter for the original 
system. This parameter can be absorbed into the defini- 
tion of fi, and does not play a role in the half-filled case. 
We simply set Vfc=o = 1 in the calculation. 

The numerical aspect of implementing the VCA has 
been studied in detaili^iiS Here we use a Matsubara fre- 
quency formalism, and the calculation is composed of 
several steps. First, diagonalize the three-site Hamilto- 
nian Eq. (66), and calculate the Green's functions G 
and n through Eqs. (71)- (73). Second, 1/M2/ and 
1/NQ.b in Eq- (69) are evaluated. Finally, we search 
for the point in the space of {e, V, £,7} where l/N&lf is 




a 
a 





I y=0.293 




/ 5*=0.875 " 




/ £2 .'=-1.964 






_ (b) - ~ - - Sj 

1 1 1 1 


1 



-0.04 -0.02 





x-x* 



0.02 0.04 



FIG. 2: The fi/ and fib in the vicinity of the metallic sta- 
tionary point at T = 0.01, U = 1.0, and V = 0.5. They 
are shifted by their stationary values fi J and fi£ , respectively, 
and magnified by f3 = 100. V , £*, and 7* denote the corre- 
sponding values at the stationary point. 



In the following, we summarize our results for the ex- 
tended Hubbard model. Since the main purpose is to 
check the theory, the study is limited to half filling and in 
the paramagnetic, charge-disordered phase. At low tem- 
peratures, a critical U c value separates the metallic and 
the Mott insulating phases. The nonlocal Coulomb repul- 
sion Vq may influence the U c value, and may even change 
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the nature of the phase transition. For Vq sufficiently 
large, a charge-ordering phase becomes more stable. The 
description of such ordered phase needs the symmetry- 
broken version of EVCA. Our discussion below resides in 
the metallic and insulating phase separately and leave 
the detailed study of the Mott transition and charge- 
ordering transition for later work. 

We set W = 1.0 as the energy unit. All the data shown 
here are obtained for a temperature T = 0.01 < T c , where 
T c w 0.014 is the critical temperature of the first-order 
Mott transition obtained for the two-site VCA. In Fig. 2, 
an example is shown for the stationary point of Slf and 
fif, at U=1.0 which is on the metal side of the Mott transi- 
tion. Since we fixed e = 0, there are three free variational 
parameters left, 7, V, and £. It is a metallic solution with 
finite V* values. The vertical scales in Fig. 2(a) and 2(b) 
shows that flf changes two orders of magnitude larger 
than f2(,, given similar variation of the parameters. This 
is a hint of why the first version of EVCA (58) works not 
so well in this case and we have to adopt the EVCA Eq. 
(62). Now we check to what extent the above three-site 
EVCA solution is close to the EDMFT solution. In the 
continuous bath limit, the EVCA Eq. (56) with semi- 
circular density of states (74) will lead to the following 
EDMFT relations: 

W 2 W 2 
A (iu n ) = —G<j {iw n ) = ~ x ° ' 



where 



$ {iv n ) - 
Y (iu n ) = 



V 2 
iSL 

4a 



II (iv n ) 



V 2 
1 



Po (c) 



A (iuj n ) — e + Ga (iuj„ ) 
Vo(e) 



-an- 1 [iv n ) - $ (iv n ) - e 



Y (iv n ) , (76) 



dc. 

de. (77) 



When the baths have only finite degrees of freedom, the 
solutions of EVCA Eq. (63) will be different from the 
continuous limit solutions, and hence Eq. (76) will not be 
satisfied anymore. Therefore the discrepancies between 
the three quantities can be taken as a criterion for the 
quality of the three-site EVCA solution. In Fig. 3(a), the 
solution of the Green's functions II (ii/ n ) at the station- 
ary point for T = 0.01, U — 1.0 is shown. The differences 
between (V^/Aa)Ii (iv n ) and $ (ii/ n ), and that between 
(V r 2 /4a)II (iv n ) and -(Vq/4)Y (iu n ) are shown in Fig. 
3(b). It is seen that the fulfillment of Eq. (76) is rather 
good for the present case, although we have only three 
variational parameters. The relative discrepancy in the 
small-frequency regime is less than 1%. This shows that 
the EVCA indeed has the similar merit of the VCA, i.e., 
it is very efficient to represent the continuous baths de- 
grees of freedom by a small cluster. For other Vq values, 
we observed that with Vo increasing, the relative discrep- 
ancies also increases. Up to the upper boundary V max of 
the parameter Vq , the relative discrepancies are less than 
10%. 



0.1 

0.08 

5 0.06 - 
F 
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— 1 — 1 — 1 — 1 — 1 — 

V 7(4a) n n - * n 




- - V 2 /(4a) (n n + aY n ) - 


1 


, 





FIG. 3: (a) The density-density Green's function at the sta- 
tionary point for T = 0.01, U = 1.0, and V = 0.5; (b) the 
discrepancies between certain quantities, which should be zero 
in the continuous limit of bath, see text, a = 1/2. 



For fixed U and T, the upper boundary V max is de- 
fined as the value of Vo above which the pole in the inte- 
gral expression of Y (iv n ) [second equation of Eqs. (77)] 
enters the range [— Vb,Vb]. This will lead to a singular 
peak in Y (iv n ) and the EDMFT equation (76) no longer 
has a solution. This problem of EDMFT was observed 
previously>i2iS In the EVCA formula, it appears in the 
second expression of Eq. (69) where the absolute value is 
used in the argument of logarithm. This does not change 
the stationary point, but can formally extend the formal- 
ism to all regime of Vq. In the regime Vq > V max however, 
we find no physical EVCA stationary point, which is con- 
sistent with the nonexistence of an EDMFT solution in 
this regime. This breakdown of the present formalism 
in the large- Vq regime is related to the instability of the 
symmetry-unbroken solutions. It is expected that in the 
regime Vq > V ma x, & symmetry-broken solution should 
be stable in the corresponding EVCA equations^ 




W74 ImG n 
ImA 



0.5 



FIG. 4: Comparison of the three quantities related to the 
electron Green's function at the stationary point for T = 0.01, 
U = 1.0, and V = 0.5. 
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In Fig. 4, we compare the three quantities 
for the electrons, W 2 /AlmG a (iuj n ), lmA(iuj n ), and 
W 2 /41mX a (iu) n ) at the same stationary point as in Fig. 
3. Since there is only one pole tv n = in the hybridiza- 
tion function A (iu> n ), the discrepancies among them en- 
large in the small-frequency limit. However, they are 
fairly small when the frequency is larger than the com- 
mon crossing point. 

In both Fig. 3 and 4, it is interesting to observe 
that there are certain common crossing points among 
the three quantities. In each subequation of Eq. (76), 
the two successive identities are not independent. One 
implies the other. Therefore the three quantities com- 
pared in Fig. 3(b) as well as in Fig. 4 always cross at 
common points. There are two frequency values of v and 
one of u>, at which the first and second equation of Eq. 
(76) are satisfied, respectively. This corresponds to the 
fact that we have two free parameters £, 7 for bosons and 
one free parameter V for fermions in the three-site unit 
cell at half filling. Therefore, the EVCA algorithm can 
also be regarded as a specific way of fixing a number of 
boson and fermion Matsubara frequencies on which the 
EDMFT equation are satisfied exactly. In the limit of 
infinite degrees of freedom of the reference system, the 
EDMFT equations are then fulfilled at all frequencies. 
This again shows that EDMFT is the continuous limit 
of one-impurity EVCA. This observation also provides a 
possibility of constructing theories parallel to the EVCA. 
One can conceive certain ways of selecting the Matsub- 
ara frequencies at which the EDMFT equation is satisfied 
exactly. The number of such frequencies represents the 
degrees of freedom of the reference system and can be 
increased in a controlled way. However, a unique crite- 
rion for the qualities of different such approaches is still 
lacking. 

For T = 0.01 and U = 1.0, the metallic symmetry- 
unbroken solution of EVCA persists up to V max « 0.72. 
With increasing interaction strength Vo, the density- 
density Green's function n (iv n ) increases uniformly on 
all energy regime, as illustrated in Fig. 5(b). This shows 
that the charge fluctuations in the metallic phase are en- 
hanced by the nonlocal interactions. At half filling, the 
total weight of n (iv n ) is given by (3 (On) — (n) 2 + 2D) = 
2/3D, where D = (n^n^) is the double occupancy. In 
Fig. 5(a), we show the double occupancy D as a func- 
tion of Vo for U — 1.0 and T = 0.01. It is an increasing 
function up to V ma x where the symmetry-unbroken so- 
lution disappears. Physically, the enhancement of the 
charge fluctuations is related to the dynamical screening 
of the local repulsion by the nonlocal interactions. This 
effect has been observed in previous EDMFT studies of 
the extended Hubbard model, 19 and may have qualitative 
influence on the metal-insulator phase transition. 34 

In the insulating phase, the charge fluctuations are 
suppressed to very much extent, and the effect of nonlo- 
cal interaction is almost negligible at low temperatures. 
Our calculation at T = 0.01, U = 3.2 gives a stationary 
point V* = 0.04, C = 1.1, and 7* = 0.012. The double 
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FIG. 5: (a) Double occupancy as a function of Vo for T = 0.01, 
U — 1.0. The vertical dashed line marks out the bound- 
ary Vmax ~ 0.72, above which the symmetry-unbroken three- 
site EVCA solution doesnot exist, (b) The density- density 
Green's function H(ii/ n ) for several Vo values. From top to 
bottom, Vb = 0.7,0.5,0.3, and 0.0. 



occupancy at this stationary point is D = 3.6 x 10~ 4 . 
It is expected that the finite V* here is due to the finite 
temperature that we use. At exactly zero temperature, 
the two-site VCA gives zero double occupancy and van- 
ishing n (iv n ). This will lead to vanishing of the Weiss 
field Hq 1 according to the EVCA equations Eq. (76). 
Therefore, we expect that for the three-site EVCA, the 
finite Vo does not play a role in the zero-temperature in- 
sulating phase, and the results are the same as in the 
pure Hubbard model. Of course, this is an artifact of 
the three-site EVCA. It is known that for the Hubbard 
model, the double occupancy D is finite even in the insu- 
lating solution at zero temperature. More detailed study 
on the effect of nonlocal interaction on the metallic and 
insulating phases and Mott metal-insulator transition is 
to be carried out systematically in the EVCA approach. 

III. CORRELATED HOPPING 

In this section, we construct EVCA for the Hubbard 
model with correlated hopping. The correlated hopping 
term is expected to play an important role in the phe- 
nomenon of itinerant ferromagnetism4£ We consider the 
following Hamiltonian, 21 

2 

ff = -EE t™ n ciP™c ja P? a +H loc , (78) 

ija m,n=l 

where P^ = n iw and P} a = 1 - n iw , and H toc = 
U Yli n it n il i s the local Hubbard interaction. In Eq. 
(78), the most general form of correlated hopping is de- 
scribed by the hopping parameter t" m . For this model, 
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the DMFTCH has already been developed^ 20 ' 21 where the 
lattice model is mapped into one impurity embedded into 
several baths. Each bath is coupled to a different com- 
ponent of the impurity degrees of freedom. The EVCA 
formula to be established in the following is an extension 
of DMFTCH to larger cluster size as well as to discrete 
bath degrees of freedom. 



A. Generalized Luttinger-Ward functional 

Similar to the construction for the nonlocal density- 
density interaction in Sec. II, we first generalize the 
Luttinger-Ward functional for the Hubbard model with 
correlated hopping in this subsection. For convenience, 
we define the projected creation and annihilation opera- 
tor as 



dt = c 1 P" 
d — r P m 



(79) 



The statistical action with a source field Gq™™ is intro- 
duced as 



S c*,c, Go 
r P rP 



J drj dr'££ 

Jo Jo ..„._ m „ 



rP 



ija mn 
mn 



- G 



(r-r') 



I dTHi oc (r) . 
Jo 



(80) 



Equation (80) becomes the corresponding action of the 
original system when Gq = Gq, and 



(Go^ir-r') 



— - n j da - t% 



5 (t - t') 



(81) 



Here, Gq 1 should be regarded as the inverse of Go in 
the entire space expanded by time, space, spin, and the 
m = 1,2 coordinates. The partition function and the 
grand potential are 



Go 



= / \{Vc* a Vc la e-^'^\ (82) 



Go 



InS 



Go 



(83) 



The Green's function G is now defined in a similar way 
as in Eq. (13), but based on the operators di m and d\ , 



GTAr-r') 
-(T T d ma (r)dt nQ (r'))s 

\{Vc^Vc ta d lma (r) d* na (r') ^ . 



Go 



(84) 



From the above definitions, it is easy to establish that 



(r - r') = -13- 





G 




6 


Gq 1 


nm 

(t> - T) 



(85) 



which will be abbreviated as 



G 



Go 



-13- 



Gq 



(86) 



The G defined above is not the usual electron Green's 
function. If we denote the one-particle Green's function 
as 

L ija (r - r') = -{T T c la (r) c] a (r'))g, (87) 
we have the following relation between L and G: 

L ija (T-T')=J2G%Z(T-T'). (88) 
mn 

It is noted that similar with the definition of {G, n} in 
Sec. II, Go is not the value of G at Hi oc = 0. This is 
because in the action Eq. (80), Go is coupled with the 
occupation-projected operators. Even if Hi oc = 0, the 
action Eq. (80) is still a complicated many-particle prob- 
lem. 

For the next step, a Legendre transformation is applied 
to f2[Go] in a similar way as in the previous section. We 
introduce an intermediate functional F, 



f g = n 
= U 

and it follows that 
SF 



Go 
"Go 

G 



-Tr 



so. 



Go 



6G„ 



G 1 



^Tr (GG 1 



SG 



G 



(89) 



(90) 



Here, the trace should include the degree of freedom of 
m = 1,2. The generalized self-energy £ and Luttinger- 
Ward functional $[G] are defined in a similar fashion, 



G n 



G-\ 



(91) 
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$ G — f3F G -TrlnG. 



(92) 



Note that the parameter a introduced in the previous sec- 
tion is now selected as a = — 1, to be consistent with the 
pure fermion field components of di ma . We get a formula 
formally same as the one of the original Luttinger-Ward 
functional, 



5^> 



G 



SG 



(93) 



and 



0SI 



G 



= $ 



G 



+ TrlnG -Tr (GG^. (94) 



The difference here is that the G, Go and £ are matrices 
in the expanded space including index m = 1,2. The 
Luttinger-Ward expression of the functional of fl is ob- 
tained immediately by simply replacing the Go with its 
physical value G , and we get 



LW 



G 



= $ 



G 



TrlnG-Tr (GG^ 1 ) , (95) 



with the following stationary properties 



sn 



LW 



G 



8G 



IG=G 



Qlw G = G 



= o, 



fi. 



(96) 



B. EVCA for correlated hopping 



In this subsection, we continue to construct the EVCA 
equation for the correlated hopping system. Doing Leg- 
cndre transformation for $[G], we get 



A 



Tr 



G 



SG 



G 



Tr 



(EG), 



and 



= -G 



(97) 



(98) 



The grand potential f2 as a functional of the generalized 
sclf-cncrgy E reads 



= A 



TrlnlGn 1 -!! 



(99) 



The grand potential functional having the desired sta- 
tionary properties is defined as 



(3Qevca 



= A 



Tr In 



(V 



(100) 



It is obtained by replacing G 1 in Eq. (99) with its phys- 
ical value Gq 1 in Eq. (81). Using Eq. (99) to eliminate 
A in Eq. (100), we get 



&EVCA £ 

1 



+ ^Trln (Go 1 



^-^Trlr^Go 1 



(101) 



Now we introduce the reference system. As usual, the 
reference system is defined on a clusterized lattice with 
the same Hi oc part as the original Hamiltonian (78). The 
action reads 



S re f c , c, Go< 
rP r 



= I ^'EE 

ija ran 
/ \ ran 

(Goc 1 ).. (r-r') 

f dr Hi oc (r) , 
Jo 



io Jo 



djria (j~ ) 



(102) 



where the Weiss field G } is block diagonal in the spa- 
tial coordinate. The generalized self-energy functional 
&evca for the reference system can be obtained along 
the same line, except that now the generalized self-energy 
£ should be replaced with £ c . By varying the parameters 
of the reference system (102), £ c can be varied through 
out its space. Here we have made use of the fact that 
the same functional form (Ievca [£] applies both for the 
original system and for the reference system. We then 
obtain 



n 



EVCA 



n 



^Trln (Go" 1 



(103) 



Using the same principle of constructing EVCA as be- 
fore, We establish the EVCA equation for a continuous 
reference system, 



sn 



EVCA 



5E C 



= 0. 



(104) 



For a discrete reference system where the Weiss field Gq c 
is parametrized by finite number of parameters t, the 
EVCA equation reads 



SQevca 



E c (t) 



St 



o, 



(105) 
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and S a pp — S c (t — 
Eq.(104) and (105), 
EVCA, 



t ap p)- Introducing Eq. (103) into 
we get the continuous version of 



G c = 
% = 



G 1 - £, 



G- 1 f~i— 1 

Oc ■ 



(106) 



Here the cluster index / is dropped. This set of equa- 
tions has the same form as the VCA equation for the 
pure Hubbard model, but with different definitions of 
G c . Here G c is the two-particle Green's function instead 
of the usual one-particle Green's function. In the case 
of one-impurity reference system, Eq. (106) reduces to 
the equation of DMFTCH proposed by Shvaika, 21 which 
is proved to be exact in infinite spatial dimension^ In 
this sense, Eq. (106) is a generalization of DMFTCH to 
cluster reference system. For the reference system with 
discrete degrees of freedom, Eq. (105) produces for each 
cluster index /, 



Tr 



G c 



5t 



(107) 



This discrete version of EVCA equation allows for more 
flexible and efficient selection of the reference system. It 
is actually a more general form, since in the continuous 
limit, the continuous EVCA equation Eq. (106) is recov- 
ered. 



IV. TRANSLATION-INVARIANT VCA AND 
EVCA 



The VCA (Ref. 24) and our construction of EVCA 
in Sees. II and III all make use of clusters with open 
boundary conditions. They can recover the CDMFT and 
EDMFT (or DMFTCH), respectively, in the case of con- 
tinuous bath degrees of freedom. However, since these 
theories give out self-energies which are not translation 
invariant, some kind of estimator should be used to calcu- 
late the lattice self-energies from the cluster ones. Here, 
we construct a different version of VCA such that it has a 
momentum-space representation similar to DCA. In the 
continuous limit, it recovers the formalism of DCA. Simi- 
larly, the EVCA formalism developed in previous sections 
can also be adapted to using periodic boundary condi- 
tions for the clusters of the reference system. In the con- 
tinuous limit, it leads to an extension of DCA formula to 
treat nonlocal interactions (or correlated hopping). 

Our discussion in the following is for a _D-dimensional 
hypercubic lattice. Other lattice structure can be dis- 
cussed similarly. The number of sites and site interval for 
the ith dimension are denoted as N and cii , respectively. 
We divide the whole lattice into geometrically identical 
clusters of size (Lia%, £2^2, Ad<zd). The number of 
total clusters is N/L = Yii=i N/Li. The site vector of a 
cluster is denoted as 1 — (liai,l 2 a2, ..^IdId), with U = 
0, Li — 1. We define the super-lattice site vector as r = 
(riLiai,r 2 L 2 a 2 , ...,r D L D a D ), with r, ; = 0, Ni/Li - 1. 
Any lattice vector n = (n\<X\, n^a^ , 7^.00.0) with n% — 
0,1,..., N — 1 is uniquely expressed as the summation 
of a super lattice site vector and a cluster site vector, 
n = 1 + r. In the first Brillouin zone of the three lattices, 
the original lattice, the cluster lattice, and the super lat- 
tice, their corresponding inverse vectors are denoted as 
k, K c , and K, respectively. We have 



In this section, we put forward the VCA and EVCA 
equations that are based on a cluster reference system 
with periodic boundary conditions. When using a ref- 
erence system with continuous degrees of freedom, the 
VCA will recover the DCA of Jarrell et ai, and the EVCA 
will become an extension of the DCA including nonlocal 
interactions. 

The DCA and CDMFT are developed independently 
to extend the original single-impurity DMFT formalism 
to finite clusters. The CDMFT is designed on clusters 
with open boundary conditions, while the DCA adopts 
periodic boundary conditions for the cluster. Both for- 
malisms recover the DMFT in the one-impurity case, and 
become exact in the limit of large cluster size. However, 
with increasing cluster size, the two methods have ap- 
parently different convergence properties. Recently, de- 
tailed comparisons between the performance of DCA and 
CDMFT have been madei^ It was pointed out that 
the real-space formula of DCA can be obtained by re- 
placing the bare electron dispersion e(k) in the CDMFT 
equations with a different one e CJ,c (K, K c ). This means 
that the formula of CDMFT can be changed into that 
of DCA by simply replacing the bare Gqq DMFT in that 
theory with the corresponding Gq DCA , and vice versa. 



k = 



/ 27rmi 
V N iai ' 



2irmD \ 
N D a D ) 
2np D 



K = 



L\a\ Ldcld 
{ 2-rrqi 2nq D 
I AW'"' N D a D 



(mi = 0,1,. ..,Ni-l), 

(Pi = 0,1, ...,Li - 1) , 

(g i = J l,...,JVi/Zi-l). 

(108) 



Similarly, any lattice momentum k is uniquely expressed 
by k — K c + K. The orthogonal and complete relations 
read 



_ \ "* ik-(m-n 2 ) _ r 
k 

l^ ei K c .(l 1 -I 2) = 5lil2) 



K, 



N 



E iK-(n-r 2 ) _ v 
e — °ri,r 2 j 



(109) 
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and 



/. We have 



in-(ki-k 2 ) 



n 



]± V p i«-(Ki-K 2 



'Ki,K 2 - 



(110) 



Here, we consider translation-invariant t ni „ 2 and assume 
that it has a periodic boundary condition on the D- 
dimcnsional lattice. Also, only translation-invariant solu- 
tions of G c and E c are considered. For lattice symmetry- 
broken solutions, one should take care of the sublattice 
indices. The VCA equation is formally the same as the 
EVCA equation Eq. (106) for the correlated hopping, 



G, 



(Ill) 



Here the Green's function is the usual one-particle 
Green's function, and Gq 1 in the frequency axis is 



Gon in2 ( iuJ n) — {lUn + /i) "nin 2 ~^ ^iiii 



(112) 



In the spatial coordinate, G c and E c are block diagonal 
and G^ 1 is not since the hopping matrix of the origi- 
nal Hamiltonian t niIl2 is not limited within the cluster. 
We write [Gq ] ni „ 2 as [Gq ]\^ 2 2 , where the upper index 
denotes the superlattice site vector, and the lower one de- 
notes the cluster site vector. The Fourier transformation 
to the super lattice index gives 

[G^] hh (K) = Y, i°o 1 ]I£ a e- <K ' (ri - ra) 

ri— r 2 

= {icj n + n)S h i 2 -e hh (K). (113) 
The matrix e (K) is defined as 

i lh (K) = j J2 * (K + K c ) e *(K+K.).(i 1 -i 2 ) j (n4) 



k, 



where e (k = K + K c ) = £ ni _„ 2 ■■».>•-■ 



e -ik(ni-n 2 ) j g 



the dispersion of free electrons. 

The same Fourier transformation on G c and S c give 
out 

G clll2 (K)= GXi 2 e " iK ' (ri " r2)=5 ci 1 i 2 . 

ri-r 2 

Y , , (K\ — vrir 2 -iK-(ri-r 2 ) _ yl 

ri— r 2 

(115) 

Here and are matrices defined on one of the clus- 
ters, which is actually independent of the cluster index 



(Go" 1 - J 
= ^E[ G o _1 (K)-S c (K) 

K 



-i 
iii 2 



iii 2 



(116) 



Finally the VCA equation Eq. (Ill) is simplified to 



G c (i^n) = jjj E [ iuJn + M - e (K) - Sc (iw„) 

K 

(iw„) = G^ 1 (iw„) - G^ - 1 (iw„) . 



(117) 



Note that K summation is only carried out in the reduced 
Brillouin zone of the super lattice. Equation (117) is 
exactly the CDMFT formula presented in Rcf. 13. The 
real-space formulation of DCA in Ref. 13 shows that 
when e(K) in Eq. (117) is replaced by a e cyc {K), Eq. 
(117) becomes the DCA equation. e cyc (K) is defined as 



CT 2 (K) = y£ e ( K + K C ) e ^ (l1 - 12 



(118) 



The periodic boundary condition of makes it possible 
to do further Fourier transformation on Eq. (117) over 
the cluster site index, and we obtain the original DCA 
equation established on the cluster K c points: 

L 



G c (K c ) = - [™n + M - e (K + K c ) - S c (K c ) 



S c (K c ) = G 0c (KcJ^-Gc (Kc)" 1 



(119) 



DCA can produce a translation-invariant cluster self- 
energy, which can be directly identified as the lattice 
self-energy, in contrast to CDMFT. 

We are interested to see what e cyc looks like in real 
space. It is found that the corresponding real-space hop- 
ping matrix element of e cyc is 

[t cv TJ 2 2 =~ J2 e(K + K c ) e iK -( ri - r2 )+ iK ^- 12 ), 



K,K C 



(120) 

The full Fourier transformation of [t cv °Y\li 2 gives an elec- 
tron dispersion 

e cyc (K, K c ) = - ]T ]T e (K + K c ) e <^-K c -K).(i 1 -i 2 )_ 

(121) 

The real-space formula of DCA can be expressed in our 
notation as 



Gc 



S c — G 0c G c 1 , 



(122) 
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where 

[G , 0DCilnin 2 = («*>„ + M) <*nm a + dn 2 • ( 12 3) 

This can be easily verified since the Fourier transforma- 
tion to Eq. (122) over the superlattice site vector will 
lead to the real-space formulation of DCA. Comparison 
ofEq. (Ill) and Eq. (122) shows that DCA and CDMFT 
can be described by a unified equation with different bare 
electron dispersion. One is e(k) for CDMFT; the other 
is e cyc (K, K c ) for DCA. The corresponding bare Green's 
function entering the theory is Go and Gqdca, respec- 
tively. Because e (k = K c ) = e cyc (K, K c ), when clus- 
ter size L increases, the number of K c points increases. 
e cyc (K, K c ) will tend to e (k) uniformly in the k space. 
This shows the equivalence of CDMFT and DCA in the 
large-cluster limit. 

The merit of the real-space formula of DCA is that it 
not only has a translation invariant form as in the k-space 
formula, but also allows for the study of the phases with 
reduced translation invariance or system without trans- 
lation invariance, such as an antiferromagnetic phase 
or disordered system. With the above discussion, it is 
straightforward to put down the VCA equations for a 
general reference system while still keeping a translation 
invariant form. This is done simply by using Gq^ca m 
the original VCA equations. We get 



EVCA with translation invariance, 



Tr 



G c 



G, 



-l 

QDCA 



SE c (t) 



St 



(124) 



This is the DCA equation adapted for discrete refer- 
ence system. In the continuous bath limit, it recovers 
the original DCA equation. 16 DCA has already been im- 
plemented with various cluster solvers ^^J^ With Eq. 
(124), it is now possible to study the discrete version of 
DCA on reference systems with smaller cluster Hilbert 
space using numerical methods such as exact diagonal- 
ization. And the flexibility and efficiency of the VCA 
is thus combined with the merit of translation-invariant 
self-energy of DCA. 

For the EVCA, a translation invariant form can also be 
obtained along the same line. For the case of the density- 
density interaction, we use the modified hopping matrix 
*n^n 2 Eq. (120), and define the interacting parameters 

ycyc ag 

rycycirjra _ _ J_ V fK + K e iK ' (n - r 2)+iK c (li -1 2 ) 

L J lila JJ / j V c > 

K,K C 

(125) 

Replacing G 1 and II 1 in Eqs. (59) and (60) with the 
corresponding DCA version Eq. (123) and 



[nooCjJnma — Kifn 2 > 



(126) 



we conveniently get the general real space formulas for 
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(pQDCA ^cj 
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ll 0DCA 1 < 



1 SE, 
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St 
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St 



= 0, 
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ODCA 
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sv 



-- 0. 

(128) 



Another version of EVCA decouples the fermion and bo- 
son contributions to the grand potential. Its correspond- 
ing translation invariant form reads 



Tr 



Tr 



G r - 



n, 



ODCA 



-E, 



6E, 



(IV) 



n 



ODCA 



St 

sf c (t, V 
sv 



= 0, 



0. 



(129) 



In the continuous bath limit, the above EVCA theories 
become the extension of DCA to include nonlocal inter- 
actions, and can be denoted as EDMFT+DCA. 

Similarly, the EVCA for correlated hopping can also be 
implemented in a translation invariant way. We simply 
replace the bare Green's function in Eq. (107) with the 
DCA version, and the obtained EVCA equation reads 



Tr 



G c 



G, 



-l 

ODCA 



5E C (t) 



St 



0. 



(130) 



where Gq^ ca (iuj n ) can be described by Eqs. (120) and 
(123), but with the t cyc and e (k) being 2x2 matrices in 
the space of single occupation and nonoccupation. 



V. DISCUSSION 

In this section, several issues about EVCA are dis- 
cussed. We first discuss the EVCA equations for the 
ordered phase. Then, the EVCA is constructed after car- 
rying out a Hubbard-Stratonovich transformation on the 
original system. Comparison with previous EVCA equa- 
tions shows that the Hubbard-Stratonovich transforma- 
tion plays a nontrivial role here. Using the same method, 
the EVCA for the Heisenberg model, a correlated spin 
system is obtained straightforwardly. The problems with 
EVCA for bosonic systems are discussed. Finally, we 
point out that our construction applies equally well to 
classical models, such as the Ising model. 
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A. EVCA for ordered phase 



where 



When the system is in the ordered phase induced by 
the corresponding nonlocal interactions, the implemen- 
tation of EVCA need to be modified according to broken 
symmetry. This issue has been discussed thouroughly 
in the context of EDMFTi 19 i 35 i 38 Here, a similar strat- 
egy can be taken to handle the ordered phase in EVCA. 
EVCA equations developed in Sees. II and III, such as 
Eqs. (59), (60), and (63), are still valid. The Hartree 
contribution should be singled out by using the normal 
ordered operator in the original Hamiltonian. However, 
care should be taken when we write down the reference 
Hamiltonian and evaluate the trace terms in those equa- 
tions. In the case of a bipartite charge-ordered phase, 
for example, the Weiss fields appearing in the action of 
the reference system should obey the two-sublattice sym- 
metry of the ordered phase. The trace terms in Eqs. 
(59), (60), and (63) should also be evaluated on the two 
sublattices. These will lead to different final expressions 
for the Qevca functional and EVCA equations from the 
symmetry-unbroken case. Detailed expression will be dis- 
cussed elsewhere in combination with model studies. 



B. The role of Hubbard-Stratonovich 
transformation 



Up to now, the EVCA equations discussed in this 
paper are constructed directly for the original electron 
Hamiltonian. One can also first carry out a Hubbard- 
Stratonovich (HS) transformation on the nonlocal in- 
teraction terms in the Hamiltonian, and construct the 
EVCA equations for the resulting effective electron- 
phonon Hamiltonian. This latter approach is widely used 
in EDMFT studies since it facilitates the employment of 
the Monte Carlo algorithm. It has been claimed that the 
two approaches are equivalent to each other Here tak- 
ing the extended Hubbard model Eq. (3) as an example, 
we study the role of HS transformation. It is found that 
in general the two approaches are nonequivalent. For 
Hamiltonian (3), the direct construction starts from the 
action (11), 



Sa[c*,c, G ,n ] 

f dr f dr' cl (r) [~G^] ija (r - r') c, ]a (r') 
Jo Jo ija 

f dr f dr' £ : n t (r) : [-H^],, (r - r') : n 3 (r') 
Jo Jo 



KcA (r) = H loc (r) + 2 £ n t (r) [-H^ 1 ] . . (n,) 

+ E n o*>*><^>- ( 132 ) 



The corresponding EVCA equations obtained for the con- 
tinuous reference system read [Eq. (55)] 



G r = 



G 1 - E ( 



U r = -a 



Gqc G c 1 r c 



(133) 



where Gq 1 and 1 are given by Eqs. (5) and (6), re- 
spectively. 

At this stage, one could apply the HS transformation 
to Eq. (131) to introduce phonon fields. This is not a 
change to the EVCA equations but a possible way to 
solve the reference system. However, there is another 
essentially different approach to construct EVCA itself. 
One can apply the HS transformation to the original 
electron action Eq. (4), obtain an equivalent electron- 
phonon action, and construct the EVCA for the nonlocal 
phonon-phonon interaction term. In this approach, after 
the Hartree contribution to the phonon term is singled 
out, we introduce the source fields for the electron and 
phonon operators as before, and obtain the action 

S B [c*,c, G , Qo] 

= f dr f dr 1 £ c* a (r) [-G^% a (r - r') c ]a (r') 
Jo Jo 



+ i 



/ dr / dT'J2^i(r):[ 
Jo Jo tj 

dr^&c^ (r)c la (t) 

f KoB (r)dr, 
Jo 



io m 



(r-r'):^(r') 



(134) 



where 



KcB (r) 



H loc (r) + 2^0, (r) [-Qo% (4>j) 
£Qol;<&><^> (135) 



and 



^o 1 



= S Do. 



(136) 



t@ i , s The actions Sa and Sg are equivalent in the sense that 

J Q locA T ' ^ ' for Qq^ = 1/4Hq c , they give the same thermodynamical 
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average values for the same operator. G is defined as in 
Eq. (13). We define the density-density Green's function, 
phonon Green's function, and phonon self-energy as 

Hsij (r - t') = (T T [m (t) - (m)} [rij (r') - (r^)])g B , 
Qij (r - r') = (T T [fa (r) - (fa)} (r') - (^)]) §B , 



(137) 



a = 1/2 since it leads to the exact solution Q — —l/2Qo, 
G = Go in the limit of the nonelectron-phonon interac- 
tion and f/ = 0. However, since the electron-phonon 
interaction is a fixed finite value in the reference system, 
this limit is actually not reachable in the EVCA calcula- 
tions. A clear criterion for the selection of the a value is 
therefore still lacking in the present theory. 



For a system defined by the action Sb, there exist the 
following equations connecting the electron and phonon 
quantities^ 



J 3h 



(138) 



Following the previous procedure, the EVCA equations 
for action Sb are readily obtained as 



G r = 



G - £< 



G, 



Qc 



T, 



G~\ T c = Q^ + aQ-\ (139) 



It is seen that the two set of equations Eqs. (131) (133) 
and Eqs. (134) (139) are similar, but the operators ap- 
pearing in the corresponding reference actions n, and fa 
have different nature. Besides Eq. (136), the only di- 
rect connection between the two sets of equations occurs 
when = 1 /4n 0c is satisfied in the self-consistent cal- 
culations. However, even at this special point where the 
two actions are equivalent, the two set of EVCA equa- 
tions are different. In such case, considering = II 
and Eqs. (136) (138), we reexpress the phonon equation 
Q c = — a|(<3o 1 ~ T c )~ 1 \ by the density-density Green's 
function II defined by Sa as 



(fioc-r- 1 ) '= (no-r 



(140) 



This equation is different from the corresponding one in 
Eq. (133). In general, Qq c — l/4±I 0c is not guaranteed in 
the EVCA self-consistent calculations, and we have even 
less connection between the two approaches. Therefore 
our conclusion is that in general there is no reason for 
the two set of EVCA equations to give out the same re- 
sults. This is true also for discrete EVCA formalism or 
for EDMFT, the single-impurity limit of the continuous 
EVCA. We should therefore take into account this differ- 
ence when we compare the results from different EDMFT 
or EVCA calculations. 

It is noted that in the HS transformation approach, 
besides the electron field we deal with real phonon fields. 
This seems to provide a justification of the selection of 



C. EVCA for correlated spin systems and boson 
systems 

Our constructions of EVCA in Sees. II and III have 
the interesting advantage of universality. It is free to 
introduce different source fields for each pair of nonlo- 
cally coupled operators, no matter it is a pair of one- 
particle operators such as c, (for the usual VCA), or a 
many-particle one such as rii (for density-density interac- 
tion) or Ci a nia (for correlated hopping). Once the source 
fields are introduced for these operators, the subsequent 
derivation procedures are the same, and similar EVCA 
equations are obtained. Except in the definitions of the 
Green's functions, the nature of the nonlocally coupled 
operators does not play a role in the derivation. There- 
fore, it is rather straightforward to apply these deriva- 
tions to systems other than correlated fermions. At the 
same time, the Luttinger-Ward functional that is most 
frequently utilized in correlated fermion systems, can now 
be generalized for correlated spin systems and boson sys- 
tems. This also holds for the EDMFT. For example, 
the EDMFT equations for the Kondo lattice model with 
RKKY interactions can reduce to that for the Heisen- 
berg model by simply discarding the electron degrees of 
freedom. As shown in the following, the EDMFT for 
Heisenberg model will be recovered by our EVCA in the 
continuous bath limit. We also discuss the EVCA for 
correlated boson systems. 

We start from the Heisenberg model 



H 



(ij) 



(141) 



where (ij) denotes summation over pairs of sites. We 
introduce the source field IT^ 1 to each component of the 
spin operators. The resulting action reads 



s,n 





= s h 



-no" 1 



(r-r')S](r'). 

(142) 



Here, Sb is the contribution from Berry's phase. The 
original system is realized when 

' (r-r') = [n -^.(r-r') 



J 13 5{t-t'). (143) 
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The partition function Z and the free energy F are 
Z = jUvSie- 8 ^}, 



F 



InZ. 



(144) 



We introduce the Green's functions 



IT7. (T-T>) = lj HVS, [SV (r) S] (/) e^'H 



and the self-energy 



r = m 



ah" 



(145) 



(146) 



Here, r\ = x,y,z are the spin components, a = 1/2 
is usually used. The latter equation is regarded as the 
matrix equation in the lattice, time, and spin space. The 
same procedure as in Sec. II will lead to the following 
form of the free energy functional: 



F, 



LW 



n 



n 



— Tr lnll — -Tr 

P P 



Ho n 



(147) 



Here, the trace is carried out in the lattice coordinate, 
time, and spin space. F^ is stationary at the physical 
value of the Green's function, and up to a constant, its 
value at the stationary point is the physical value of the 
free energy. The derivative of the generalized Luttinger- 
Ward functional ^[n] with respect to n gives the self- 
energy. The frcc-encrgy functional necessary for the con- 
struction of EVCA equation reads 



Fevca 



F 



r 

1 

0' 



-Tr In H 



1, 



-Tr In H 



(148) 



From this functional, the EVCA equations for the Heisen- 
berg model are obtained as 



n c 

Tr : 



n 



n 



-i _ 



all- 1 



r, 



(149) 



Equation (149) is formally the same as the general VCA 
equation for electrons, but with the reference system de- 
fined as 



Sref 
P 



s,n 0( 

■p 



ija 



Sr(r)\-U^] a (r-r')Sf(r'). (150) 



For a discrete reference system, n^ 1 can be parametrized 
by some parameters J, and the first equation in Eq. (149) 
is replaced with 



Tr 



n, 



IK 



sr c (t) 
si 



o. 



(151) 



With this equation, it is possible to study the Heisenberg 
model with discrete reference systems. 

The EDMFT has been used to study the coupled 
fermion-boson system^ where the bosons are phonons. 
Recently the correlated boson systems is studied inten- 
sively in the context of cold boson atoms in the opti- 
cal lattice. The bosonic Hubbard model^ is one of the 
simplest models to describe the boson atoms which have 
particle number conservation. 



loc 



(152) 



where Hi oc = UJ2i( a t a i ~ l) a l°i describes the local re- 
pulsion between bosons. The ideal EVCA for this model 
should be able to describe both the Bose-Einstein con- 
sensation phase and the Mott insulator phase. For this 
problem, it is important to know that what kind of Weiss 
fields are necessary for the description of the physics. 
From the usual splitting of the boson field into sum of 
condensate wave function and the quasiparticle field, it 
is plausible that besides the normal boson Green's func- 
tion, the new Weiss fields should include the condensate 
fields as well as the anomalous Green's functions. 46 The 
corresponding self-energy functional and EVCA equation 
can be derived in a similar way as for fermions. Detailed 
work in this direction is in progress. It would be interest- 
ing to explore the EVCA for correlated bosons and make 
comparison with the results from other approaches. 



D. EVCA for classical systems 

In this part, we derive the EVCA equations for a clas- 
sical system, the ferromagnetic Ising model. The purpose 
is to show that EVCA applied on larger clusters may be 
used as a useful tool for the study of classical systems. 
Here we consider the classical Ising model 



(153) 



Sf is the z component of spin 1/2. The summation is 
over pairs, and we assume Jy > 0. Because all the op- 
erators present in this model commute with each other, 
there is no quantum fluctuations in this model. In this 
case, instead of introducing a time-dependent source field 
as for the Heisenberg model, it is sufficient to introduce 
static ones [n^ 1 ]^ . Taking into account of the possible 
symmetry-broken phase, the reference system with an 
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open cluster boundary condition has the action 

Sref \S Z , Iloc 



J tj (S-mj+S-rrii 

(ij) 



mirrij) 



and the partition function 



Z = Tre 



(154) 



(155) 



Following the same procedure as before, we get the fol- 
lowing EVCA equations: 



n, 



= MS?-m i )(S?-m j )) Sr .„ 



Tl^ 1 +aU~ 1 , 



—a 

(Si 



n 



Qij 



2 Jij- 



(156) 



Equations (154)-(156) form a closed set of self-consistent 
equations for the Ising model. The Weiss fields are real 
variables. 

For the simplest case of the one-site cluster, we obtain 
the following self-consistent equations: 



1 



4 

-E 

2N ^ 



-1/2.4 



l/[2/3(l/4-m 2 )] 



tanh 



P n - x + -J m 



(157) 



This set of equations is different from the usual Weiss 
mean-field equations. It has been obtained (with differ- 
ent parameter definitions) by Pankov et al. in a semi- 
classical analysis of EDMFTii It was found that this 
mean-field theory predicts a better phase transition tem- 
perature than the Weiss mean-field theory. However, the 



order of the transition is given incorrectly. Instead of a 
second-order transition, this one-site theory gives a first- 
order phase transition for dimensions d < 4. It is ex- 
pected that the approximation will be improved by using 
a larger cluster size. However, it is still an open question 
whether the correct order of the phase transition can be 
obtained in EVCA with larger clusters. 



VI. SUMMARY 

In this paper, we extend the VCA to include the nonlo- 
cal interactions. A universal way of constructing this the- 
ory is described for different systems, i.e., the Hubbard 
model with density-density interaction, the Hubbard 
model with correlated hopping, the Heisenberg model, 
and the Ising model. The construction method proposed 
here is universal, and may deserve further consideration 
in the context of the general variational principled Tak- 
ing the extended Hubbard model as an example, we show 
that the simplest three-site EVCA can already produce 
results very close to those of the EDMFT. The latter can 
only be realized in the limit of a continuous reference 
system. A number of existing theories are obtained as 
limiting cases of EVCA. We arrive at a cluster extension 
of EDMFT for a system with nonlocal interactions, and 
arrive at a cluster extension of the DMFTCH for systems 
with correlated hopping. VCA (EVCA) equations which 
are based on clusters with periodic boundary conditions 
are also proposed. This produces translation-invariant 
solutions and has the DCA (DCA extended to include 
nonlocal interaction) as the continuous limit. The high 
efficiency and flexibility of the EVCA makes it possible to 
study a large class of correlated systems with less numer- 
ical cost by adopting a relatively simple reference system. 
Some interesting issues about EVCA are discussed. 
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